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We analyze the propagation of fast-light pulses through a finite-length resonant gain medium both analytically 
and numerically. We find that intrinsic instabilities can be avoided in attaining a substantial peak advance 
with an ultra-short rather than a long or adiabatic probe. © 2008 Optical Society of America 
OCIS codes: (270.5530) Pulse propagation and solitons; (020.1670) Coherent optical effects 



Studies of electromagnctically induced transparency 
(EIT)^ have suggested ways that a strong laser field can 
dress a dielectric medium such that a weak probe pulse 
injected in the dressed medium will exhibit non-standard 
velocity of propagation.^ Many observations of so-called 
slow light'^ and fast light'' in such dressed media have 
been reported. Here the terms fast and slow refer to 
advanced or retarded propagation of a pulse peak, giv- 
ing the appearance of pulse propagation faster or slower 
(even much slower) than the group velocity in the same 
medium without the benefit of dressing. 

In the case of fast light, which is our principal con- 
cern here, an effective technique^ is the manipulation, 
by strong coupling fields, of the spectral domain be- 
tween two resonance lines by external cw fields to obtain 
very high dispersion in a very narrow spectral window 
through which a probe pulse will propagate with strongly 
altered velocity and anomalously low loss. The narrow 
spectral window can be obtained in more than one way^ 
but it is a universal feature of recent work and it severely 
limits the bandwidth of the probe pulse. Noise-related 
limitations on fast light have also been discussed.^ We 
have considered modifications that may permit narrow- 
line limitations to be overcome, and we report results of a 
theoretical study of fast light with very short rather than 
very long probes. Fast light propagation of short pulses 
appears feasible with the use of appropriately configured 
transient gain media, given sufficient attention to natural 
instabilities. 

Fast light is associated with negative group velocities 
(for reviews and accounts of early experiments see Refs. 
7,8). Before 1990 theoretical examinations of this effect 
in gain media were largely semi-quantitative or numer- 
ical (e.g., see 9, 10), but propagation studies of super- 
fluorescence by inverse scattering methods gave a num- 
ber of valuable analytic insights. ^^ Here we have used 
an approach that provides analytic expressions for fast 
probe pulses and is flexible enough to treat different sit- 
uations of potential experimental interest, including seg- 
mented and stacked dielectric media. Our solutions to 
the coupled nonlinear equations governing near-resonant 
pulse propagation are valid when probe pulses and tran- 
sit times are short enough to allow upper-level relaxation 



effects to be ignored. As we mentioned above, this is out- 
side the zone of KIT, which uses long pulses and rapid 
upper level relaxation to produce a quasi-steady state, 
but we will show that the pulse advance associated with 
fast light can still be produced. 
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Fig. 1. Two-level atom with level 1 connected to level 2 
via the Rabi frequency fJ of a laser field detuned from res- 
onance by an amount A, and sketch of significant pulse 
advance in an inverted two-level gain medium. 

To be as simple as possible, we model the dressed 
medium by a collection of inverted two-level atoms, sim- 
ilar to an approach taken by earlier workers.^' '^'^ A typ- 
ical two-level atom is shown in Fig. ^ and we consider 
media nominally subject to inhomogeneous broadening. 
The Hamiltonian of such a system in the rotating wave 
picture is given by 

H = -h^\l){2\-h^\2){l\ + h^2){2\ (1) 

where J7 = 2fiE/h is the Rabi frequency associated with 
the 1 — > 2 transition and for simplicity it is assumed to 
be real, /i is the dipole matrix element, E is the electric 
field envelope function, and A = uj2i —lu is the detuning 
of the laser below resonance. 

We assume an atomic state of the form \ip{x,t)) = 
ci{x,t)\l) + C2{x,t)\2) with ci and C2 being the corre- 
sponding probability amplitudes. Schrodinger's equation 
along with the Hamiltonian in Eq. Q then gives 
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Maxwell's equation written in the slowly varying enve- 
lope approximation (SVEA) yields an equation for the 



pulse envelope via the Rabi frequency: 
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and Tn* is the inhomo- 



where F(A) = ^e-^^^*^)'/^, aii^ ., 
geneous lifetime. The product ■KgF{0)/c is the inhomo- 
geneously broadened Beers absorption coefficient, and 
g = Nii^ui/eoh, N is the density of resonant atoms and 
Lu is the probe pulse frequency. 

The key differences to EIT theory are evident here - 
there are no empirical relaxation terms included in the 
equations above, and there is no coupling laser to a third 
level. The first is justified by our desire to work with very 
short pulses, so we are thinking of pulses shorter than up- 
per level dephasing times and lifetimes of say 10 ns. The 
second is based on the assumption that the conceptually 
simplest dressing of the medium is partial or complete 
initial inversion, which we will see to be sufficient for 
our purpose. We will see that the consequences of the 
differences are substantial. 

We focus here specifically on stimulated rather than 
spontaneous emission, so all interactions are taken as in- 
duced, assuming negligible influence from spontaneous 
emission on the time scales of interest, including super- 
fluorescence,^^ which merits a speciflc qualiflcation, men- 
tioned later. 

We solve Eqns. Q and ((SJ for a medium located be- 
tween the flnite positions xq and xi. The method we 
use, a variant of the Backlund transformation technique 
developed by Park and Shin,^"^ allows us to flnd solu- 
tions for media that have a variable density of atoms. 
We have taken the simplest special case, constant A^ in- 
side the medium and N = outside, as is relevant to a 
thin inverted medium. 

Before the probe pulse arrives, the state of the atoms 
in the medium is \tlj{x, t = — oo)) = |2). The exact SVEA 
pulse solutions are given by the following segmented ex- 
pression that describes a fully continuous pulse in time 
and space: 



n{x,t) 



fsechi(t-f) 
^sechi (t-^ 
^sechi it-s. 



X < Xf) 

Xq < x < Xi (4) 

X > Xi, 



where r is the nominal probe peak width. The group 
velocity Vg is given by 



Va 2 



^(A)- 



dA 



{l/rf 



^, (5) 



where the final limit applies when T2 ^ t. In general, 
Vg > c or Wg < 0, depending on T2*, the width r of the 
probe peak and g the atom-field coupling parameter. The 
expressions 4>q = —T^gx^j^c and 0i = T^g{xi — xo)/2c 
are given in the long Tj* limit. The corresponding proba- 
bility amplitudes for the atoms in the medium are found 



to be given by ci{x,t) = i SGch.(j){x,t) and €2(2:, t) 
— tanh0(a;, t), where 0(a;, t) = i 1 1 — -^ + 00 
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Fig. 2. Snapshots of the analytic solutions given in Eq. 
Q). The horizontal axis is x in units of ct. The shaded 
region indicates the location of the medium, whose fea- 
tures are described in the text. For an EIT comparison, 
see Fig. 15 in Ref. 16. 

Pulse and amplitude solutions for a fully coherent fl- 
nite length medium were previously obtained by An- 
dreev^** in the absence of Dopplcr broadening. ^^ An adi- 
abatic EIT-based example has been given by Milonni.^^ 
Snapshots of our solution in Eq. Q) are plotted in Fig. 
121 where one can clearly see that the pulse entering the 
medium is advanced - in the third frame an outgoing 
peak leaves the medium before the input peak enters 
it. In terms of initial peak widths this peak advance 
is slightly more than 10 peak widths. In this particular 
case, the sample was chosen to be Rb vapor, using the 
D2 line with these parameters: Tj* = 0.733 ns, g = 266 
ns~^ and r = 0.1 ns, all of which appear to be within 
experimental reach. They give an inverse Beer's length 
of a = yJiT /2gT2 /c w 8.15 cm"-'^ and a negative group 
velocity in the medium oivg/c ~ —3.27, and the medium 
was taken approximately 250 absorption lengths deep. 

The analytic pulse solutions Q) have 2tt Area, and 
thus after interacting with the medium the pulse returns 
the atoms to their excited state. This is a fully coherent 
but unstable situation, and the trigger for instability can 
conic from several mechanisms. By assuming a very short 
pulse we have eliminated the main ones, which are asso- 
ciated with the medium's relaxation processes. However, 
we cannot realistically assume infinitely long pulse wings 
(as in the analytic sech functions), and the finite front 
edge of a pulse is an instability trigger that will always 
be present even if very weak (see 17, 18 for a discussion 
of the consequences) . 

To test sensitivity to this semiclassical trigger of insta- 
bility we have generated numerical solutions to the same 
set of equations, with the same medium parameters. We 
used a short enough probe peak to justify ignoring sig- 
nificant medium relaxation, but we gave the pulse finite 



weak leading and trailing edges. The results are shown 
in the following figure. 
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Fig. 3. Snapshots during evolution of the numerical 
pulse solution to Eqns. jSJ and J^l with the same pa- 
rameters as for Fig. |21 but for a pulse initially entirely 
outside the medium with abrupt edges at ± 25 units 
from the peak. The pulse advance is 10 units forward, 
thus still well behind the front edge of the pulse. 

In the snapshots of the numerical solution one again 
sees in frame 3 an advanced peak leaving the medium be- 
fore the input peak enters, and there is a clear negative 
group velocity in the medium. The output pulse is rela- 
tively cleanly separated from background, and is essen- 
tially identical to that for the analytic solution, in both 
its extent of advance and in its Area, which is very close 
to 27r. The pulse "ringing" trails the advanced pulse is ra- 
diation arising from the stimulated coherence imparted 
to the medium by the leading edge of the pulse. The 
pulse must have Area = ±7r to satisfy the two-level Area 
Theorem^^ for a gain medium, and numerical integration 
confirms that the total Area of the pulse changes from 
27r to TT upon propagation through the inverted medium. 
Burnham-Chiao ringing^" is not visible in Fig. |2| but 
appears in later frames. 

True spontaneous effects such as supcrfluorescence are 
excluded in any treatment using classical fields but we 
have made a preliminary study of spontaneous effects on 
pulse advance by testing propagation in inverted media 
with small randomly distributed zero-average dipole mo- 
ments, as a surrogate for the vacuum fluctuations that 
"induce" spontaneous emission. These results indicate 
significant observable pulse advance (i.e., greater than 
1 full input pulse width) even in the presence of such 
fluctuations. 

In summary, we have examined a regime of semiclas- 
sical pulse propagation opposite that of EIT. To do this 
we chose pulse peak durations very short, and neglected 
medium relaxation times, so quasi-steady or adiabatic 
atomic states could not develop. A substantial peak ad- 
vance is evident in our exact short-pulse solution, which 
has 27r Area and infinitely extended tails on both sides 



of its brief intense peak. With inverted media, such ex- 
act solutions are well known to be unstable. The exact 
solution pulse served as a template against which our nu- 
merical propagation in equally inverted media, but with 
more realistic pulses with finite tails, could be compared. 
Our results indicate that a significant pulse advance (i.e., 
the fast light effect) appears reasonable to achieve, well- 
separated from the instability. In fact. Fig. 13 shows an 
equal advance in the numerical pulse as in the analytic 
result. 

It is interesting to note that since r ~ T2*/7.33, 
Doppler broadening played a limited role in this ex- 
ercise. An increased inhomogeneous bandwidth would 
allow more rapid dephasing of atomic coherence and add 
a degree of stabilization. 
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